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Abstract 



Q^l The phonon spectrum of Coulomb lattice in neutron star crusts above 

■ the neutron drip density is affected by the interaction with the ambient 

, neutron Fermi-liquid. For the values of the neutron-phonon coupling 

Cn| I constant in the range 0.1 < A < 1 an appreciable renormalization of the 

phonon spectrum occurs which can lead to a lattice instability mani- 
fested in an exponential growth of the density fluctuations. The BCS 

ON ' phonon exchange mechanism of superconductivity leads to neutron pair- 

en ' 
(N 



ing with a gap in the neutron excitation spectrum comparable to that 
due to the direct nuclear interaction. 



o 

OO 

Cn ■ 1 Introduction 

r-j 

Below the melting temperature T m ~ Z 2 e 2 /100rj ~ 10 9 — 10 10 K, where fj is 

6 



the average interion spacing, Z is the ion charge, the ionic component of plasma 
in neutron star crusts is arranged in a Coulomb lattice. The screening of the 
electrostatic potential of the lattice by electrons is ineffective since the ratio 
Ao/r, < 1, where Ad is the Debye screening radius, and, therefore, electrons 
are distributed almost uniformly. Above the neutron drip density 4 x 10 13 
g cm -3 the intervening space between the clusters is filled by a neutron fluid 
which goes over to a superfluid state below the critical temperature of the order 
T c ~ 10 9 K; (for a review of the early work see, e.g., ref. 0; recent progress 
is summarized in ref. @). Although initial studies of the collective effects in 
the crusts were focused either on the properties of the highly compressed solid 
matter or the unbound neutron fluid at subnuclear densities (with an exception 
of the equation of state where phase equilibrium conditions are imposed) this 
separation is justified only when the interaction between these components of 
the crusts is weak; this is not, however, always the case. The purpose of this 
work is to continue the discussion of the interaction between the neutrons in 
the continuum and the excitations of the crustal lattice, i.e. phonons, set up 
in an earlier work [Bfl. 




Figure 1: The phonon Dyson equation and the vertex equation for coupled system of 
fcrmions (solid lines) and phonon (wavy lines). The fcrmion Dyson equation is given by the 
block in Fig. 3. The thick lines correspond to the full propagators, while the thin lines to 
the free ones. 



2 The Phonon Spectrum 

The self-consistent coupled fermion-phonon problem is shown in terms of Feyn- 
man diagrams in Fig. 1 and 3. In this section we shall give the finite tempera- 
ture solution of the Dyson equation for phonons which are coupled to fermions 
of arbitrary relativism; in the present case the neutrons are only mildly rel- 
ativists (if at all) while the electrons are ultrarelativistic. The problem can 
be solved for an arbitrary coupling strength since the perturbation series con- 
verge rapidly with respect to another parameter - ratio of the phonon to the 
fermion energies: as a result one may consider only to the lowest order term 
in the integral equation for the vertex (Migdal theorem). The sequence in 
which the Dyson equations for the fermions and phonons should be solved is 
fixed by the fact that the phonons affect only a narrow range of energies in the 
fermion propagator (of the order the Debye frequency lvd), which implies that 
the polarization function can be evaluated with free fermion propagators: 

%uj v - (e p - pi) 2 { e p e p J 

where e{p) = y/p 2 + m 2 , a = 707 and we have kept only the positive energy 
states from the outset; (to excite an excitation at the top of the Fermi sea and 
an anti-particle at rest one needs to overcome an energy barrier of the order 
of the Fermi energy €f, however u < ud <C e^). The Matsubara frequencies 
assume discrete odd integer values io v = (2u + l)/(3 where v = 0, ±1, ±2, . . . 
and (3 is the inverse temperature; the other notations have their usual meaning. 
We find the retarded polarization function depicted in Fig. 1 by performing 
the frequency summation over the product of the fermion propagators and an 
analytical continuation to the real axis: 

u^(k,.) = Ml /-| V (£p+fc " ^ " f(€p tS ^ 1 + — } 
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Figure 2: The renormalization of the phonon spectrum u>(k) and the damping j(k) of 
the modes (light lines) as a function of the momentum transfer for several values of the 
neutron-phonon coupling constant. 

where M.k is the fermion-phonon coupling matrix element and f(x) = [1 + 
exp(/3x)] _1 is the Fermi distribution function. The scalar polarization func- 
tion above corresponds to particle-hole excitation for arbitrary degree of rel- 
ativism of the system and finite-temperatures; for applications to the neutron 
subsystem one may take the non-relativistic limit. Fig. 2 shows the real part 
(more precisely 1 + ReII(k, u)) and the imaginary part of the polarization 
function as a function of the momentum transfer for (3~ l = and a typical 
ratio c s /vp = 0.1, where c s and vp are the sound and Fermi velocities, respec- 
tively. The first quantity gives the renormalization of the phonon frequencies 
[u(k)/u>o] 2 (where c^o is the unperturbed phonon frequency) which are treated 
in the Debye model. The eigen modes of the system vanish first in the long- 
wave limit where, with increasing coupling constant, the curves cross zero. The 
disappearance of the real solutions to the dispersion relation indicates insta- 
bility of the system with respect to the density fluctuations, as can be seen by 
inspecting the proper solutions which exhibit exponentially growing amplitude 
of phonon mode. 



Fi gure 3: The Dyson equations for the normal and anomalous fermion propagators. The 
contributions from the anomalous sector are shown by the light lines. 

If this is the case, one may conclude that the starting lattice structure does 
not correspond to the true minimum of the energy of the system (in other 
words is not the 'true vacuum'). Available estimates of the neutron-phonon 
coupling constants are, however, unreliable for two reasons: first, preliminary 
considerations [|3] , based on a fit to the elastic neutron-nucleus scattering cross- 
sections, does not include the medium effects (which most likely will reduce the 
effective cross-section) and, in addition, the substantial imaginary part of the 
optical potential responsible for the absorptive processes should be included in 
the estimates. Second, the compositions of the crustal matter due to different 
authors (see refs. in 0, ||) imply coupling constants from marginal up to 
of the order unity at no/3, where no is the nuclear saturation density, and 
consequently, very different assessments about the role of the neutron-phonon 
interaction. 

3 Phonon Exchange Interaction and Neutron Pairing 

The net two-body interaction between neutrons in the continuum comprises the 
nuclear component (which is dominated by the attractive 1 Sq channel at densi- 
ties of interest) and the component due to the phonon exchange. Since the lat- 
ter interaction is attractive as well, the BCS pairing for some values of neutron- 
phonon coupling constant might be influenced by the phonon-exchange mech- 
anism. The respective gap equations, which include the retardation of the 
effective interaction, emerge as a solution to the set of diagrams in Fig. 3. We 
find (for simplicity the fi^ 1 = limit is given below) 





where z{uj) is the wave function renormalization and the effective interactions 
via phonon exchange and direct nuclear force are defined respectively 
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Here the nuclear interaction is given by the time-local interaction V(k); note 
that a pairing description based on the pion exchange model would allow for 
the effects of retardation, pion mode softening, condensation etc. The formal 
structure of equations (g) coincides with the Eliashberg equations for metallic 
superconductor in the presence of the repulsive Coulomb force between the 
electrons. For the values up = 0.7 MeV and A = 0.45 (which correspond to 
the density 0.3n , for the composition of Arponen, see J|) and W = 0, eqs. 
@ predict a value A = 0.03 MeV, which is basically less, but comparable to 
the gap found neglecting the neutron-phonon coupling. 

To conclude, if in a certain region of neutron star crusts the neutron-phonon 
coupling is not negligible (i.e. A < 1), a number of novel effects emerge, 
which are potentially important for the correct description of phenomena in 
neutron star crusts. The crucial problem is whether the calculations of the 
coupling constants beyond the simple estimates would lead to significant values 
of the neutron-phonon coupling constant. In addition one may need a better 
understanding of the sources of the discrepancies between different models 
of neutron star crust composition and a reassessment based on the progress 
achieved in the many-body theory of nuclear matter in recent years. 
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